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00 ; Abstract 

o 

Azimuthal anisotropy in heavy-ion collisions is studied by taking into account the 
ridges generated by semi-hard scattering of intermediate-momentum partons, which 
can be sensitive to the initial spatial configuration of the medium in non-central colli- 
sions. In a simple treatment of the problem where the recombination of only thermal 
partons is considered, analytical formulas can be derived that yield results in accord 
with the data on V2 for px <l-2 GeV/c. Centrality dependence is described by a ge- 
ometrical factor. Ridge phenomenology is used to determine the initial slopes of vi 
at low pt for both pion and proton. For higher px, shower partons from high-py jets 
must be included, but they are not considered here. 
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1 Introduction 



Hard scattering of partons is known to produce low-p^ hadrons by fragmentation, but it oc- 
curs so rarely that its influence on the bulk medium in heavy- ion collisions is not important. 
Semi-hard scattering, however, occurs much more frequently. Although there is no reliable 
way to calculate weak jets, their effect on \ow-px phenomenology should not be ignored. 
Recent experiments at RHIC have shown that hadrons associated with triggers at interme- 
diate px exhibit ridge structure at low px that protrudes above the background [H, El E]. 
Such structure has also been seen in auto-correlation without triggers pE]. The question that 
we address here is what aspect of the properties of the hadrons produced at low px can be 
influenced by the ridges generated by the semi-hard scattering. Since the time scale involved 
in such scattering is short, the ridges that are produced near the surface are sensitive to the 
initial configuration of the medium. We show in this paper that such semi-hard scattering 
processes can drive the elliptic flow of the medium. The mechanism represents an alter- 
native to the conventional hydrodynamical description that relies on rapid thermalization. 
Our approach does not replace hydrodynamics in general, but releases it from the burden of 
assuming its validity at early time. 

We restrict our attention in this paper to the low px region with px < 2 GeV/c. The 
hadronization process will be treated by recombination, as we have done in the region px > 2 
GeV/c [5]. In the higher px region the shower partons from high-p-r jets must be considered; 
their contribution to v 2 will be investigated in a future study. That region is not germane to 
the issue at hand of relating the ridges to elliptic flow at low px where only thermal partons 
are relevant. The ridges we consider are generated by semi-hard jets that are produced much 
more copiously than the high-momentum jets. The formation of ridges with or without 
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identifiable peaks has been studied before in the recombination model [6]. We now apply 
similar consideration to the study of azimuthal anisotropy. The relevance of recombination 
has been demonstrated by the observation of quark number scaling of v% [3,[H]- 

2 Ridge Properties 

Let us adopt the very simple description of nuclear geometry with the assumption that the 
colliding nuclei have sharp boundaries at the spherical surface of radius Ra- At impact 
parameter b the almond-shaped overlap region in the transverse plane is bounded by two 
circular arcs whose angular ranges are 

|0| < $ and \tt — (p\< $, (1) 

where <fi is the azimuthal angle measured from the centers of the two nuclei, respectively, 
and 

$ = cos- 1 (6/2i? A ). (2) 

That is the spatial configuration of the system at the initial time immediately after collision. 
The momentum of a detected hadron will have an azimuthal angle also denoted by <ft. The 
expansion of the system at later times may well be described by hydrodynamics, say, for 
t > 3 fm/c, so by the time of hadronization the medium has characteristics of flow. Our 
emphasis here is to replace early time hydrodynamics by ridge phenomenology. 

The issue concerning the dependence of the detected hadrons is about the physics that 
is sensitive to the initial configuration of the system. If one considers semi-hard scattering 
as the virtuality is lowered, there comes a region when the process is soft enough to have 
frequent occurrences but hard enough to create intermediate-pr jets at early times. If those 
jets are in the 2-3 GeV/c range, the corresponding time scale is r < 0.1 fm/c, short enough 
to be sensitive to the initial spatial configuration. 

In events triggered by intermediate-pr particles it has been found that jet structure 
consists of a prominent peak above a ridge in the At] distribution [lj. Quantitative study of 
the relative strength of the peak (denoted as Jet, or J) and the ridge (R) shows that J/R 
becomes small when the trigger momentum is low and the associated-particle momentum is 
even lower (< 2 GeV/c) [2j[3]. When the R yield dominates over the J yield, the underlying 
jet is hardly recognizable and has been referred to as the phantom jet [9], since it is a Jet-less 
jet. The ridge has been described as the recombination product of enhanced thermal partons 
due to the energy lost by the semi-hard parton traversing the medium [6|. The ridge being in 
the immediate vicinity of the jet direction implies that the semi-hard scattering takes place 
near the surface of the medium and that the thermalization of the energy lost is restricted 
to the region near the jet trajectory. 

Applying the notion of phantom jets to our present problem, we have a collection of 
those semi-hard jets on the surface of the almond-shaped region, if the virtuality of the jets 
is small enough so that many of them are produced in a nuclear collisions. At t] ~ the 
semi-hard scattering involves low-a; (< 0.03) partons, which are abundant. The direction 
of a scattered parton is random, but in any given patch of spatial region near the surface 
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the average direction of all outward partons is normal to the surface, since that is the only 
direction in the geometry of the problem. When integrated over the surface, there is a layer of 
emitters of semi-hard partons whose directions are prescribed by the geometry of the overlap 
at the initial time. The ridge of enhanced thermal partons may develop subsequently, after 
allowing some time for local thermalization to take place near the surface. But the average 
directions of the phantom jets are specified by Eq. (pQ), independent of the time it takes for 
hadrons to form. Since those phantom jets are the products of semi-hard partons directed 
outwards, the recoil partons directed inward are absorbed by the medium and contribute to 
the thermalization of the bulk that has no preferred direction of expansion in the transverse 
plane. In the present study we ignore any aspect of the longitudinal expansion and consider 
only hadrons at mid-rapidity; thus the property of the ridges in t] is not as important as their 
distribution. In summary the system under study can be described by two main types of 
thermal partons that are generated by the collision: (a) the bulk that is isotropic in 0, and 
(b) the ridges that is nonzero in the region specified by Eq. ([T]). 

In Ref. [5] it is shown how the pion distribution can be obtained from the parton 
distribution in the recombination model. The region of interest there is for px > 2 GeV/c, 
so shower partons are important in addition to the thermal partons. We now use the same 
formalism for px < 2 GeV/c and consider only the thermal partons. For the bulk medium 
the invariant distribution of the thermal partons just before hadronization is 



where, with our attention focused on only the momentum in the transverse plane at y = 0, 
qx is the parton transverse momentum and go is its energy. The value of C is of no concern 
here since it will be canceled in a ratio to be calculated. Adopting the same formalism of 
recombination described in [5], we have for the pion distribution with the pion mass being 
neglected 

riN B r 2 

which describes the observed low-pr behavior of the bulk (B) medium having the parton 
distribution given by Eq. In (j4]) there is no dependence. 

In addition to the thermal partons associated with the bulk medium described above, 
we now consider the enhanced thermal partons associated with semi-hard jets. A semi-hard 
scattering that can lead to a ridge must occur near the surface because otherwise the scattered 
parton would be absorbed in the interior of the medium and never emerge from the surface 
to generate any observable effect distinguishable from the isotropic background. The average 
direction of the scattered parton is normal to the initial surface, so let us focus on such a 
representative semi-hard parton. Being massless it travels at c and can reach the expanding 
surface before that surface deviates too much from the initial configuration. However, that 
is not when hadronization occurs. What happens to the semi-hard parton after it leaves the 
medium need not concern us here, since our present focus is on hadrons with px "z 2 GeV/c. 
The emergent semi-hard parton will generate shower partons that lead to hadrons at px > 2 
GeV/c. Of central importance to us is that in traversing the medium the semi-hard parton 
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loses energy which leads to local enhancement of the thermal energy of the medium. Thermal 
equilibrium may or may not be established in less than 1 fm/c, but energy dissipated by 
the semi-hard parton that passes by has its effect in the formation of ridge downstream. 
The medium must expand whether or not it can be described reliably by hydrodynamics 
starting at To = 0.6 fm/c. By the time the density is low enough for hadronization to 
take place the geometry of the boundary of the medium may be very different from that 
of the initial configuration. Nevertheless, the region of the enhanced thermal partons is 
near the flow direction associated with the initial passage of the semi-hard parton, which is 
restricted by (1). We know from ridge phenomenology that the ridge particles are strongly 
correlated in with the trigger particles [2]. Although we have no trigger in the present 
problem, the semi-hard parton that generates a ridge specifies the direction to which the 
ridge particles are correlated. Moreover, we expect a large number of semi-hard scattering 
to produce a continuum of ridges. Particles in that continuum are therefore in the range 
of 4> corresponding to the angular window in which the semi-hard partons are emitted with 
some round-off at the edges. The emission of semi-hard partons is not to be confused with 
the formation of soft hadrons in the ridges, which occurs much later; however, the ranges of 
<p of the two are strongly correlated. 

There are many aspects in this non-perturbative problem that cannot be formulated 
rigorously. In the first place the virtuality of hard scattering is not narrowly restricted by 
any observable. Our main criteria are that it should be soft enough so that the yield of 
scattered partons is high, but hard enough so that the time of scattering is short. The 
depth from the initial surface where the semi-hard scattering occurs can also vary over some 
range. The amount of energy loss depends on both the virtuality and depth, as well as 
on the density of the medium near the surface, which is not uniformly distributed in 0. 
The non-uniformity is particularly severe near the tips of the almond-shaped overlap region. 
Thus the properties of the ridges formed at <fi near $ may be quite different from that near 
the reaction plane, which is where the most important contribution to the second harmonics 
in cf) arises. The many complications that cannot be described rigorously force us to make 
simplifying assumptions that make possible our treatment without losing the essence of the 
problem. To that end, we make a sharp cutoff in at $ for the ridge particles in the same 
spirit as in approximating the Woods-Saxon distribution by a sharp edge. We further use 
just one parameter X" to summarize the uncalculable dynamical process of ridge formation, 
and write 



This simple description contains two main features: (a) it is a consequence of the semi-hard 
jets near the surface at early time, and thus has the <fi dependence characterized by the 
pre-equilibrium partons; (b) the exponential qx dependence follows from the assumption 
of thermalization at later time, as in ([3]). The sharp edges in (jSj) can be made smoother 
and the region in between may depart slightly from uniformity depending on density vari- 
ation, but the integrated effect of the anisotropy in should not depend sensitively on the 
approximations made, which have the virtue of rendering the final result in closed form. 

Both T and T' may be characterized by hydrodynamical flow and are properties of the 
pre-hadronization partons. We do not treat the problem between the initial and final states. 



dN q B+R 





(5) 
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We take the values of T and X" from the data to be discussed below. Applying recombination 
to Eq. (jSj) as before, we obtain the pion distribution for both B and R to be, for \<p\ < $ 
and 1 7T - <p\ < $, 

B M + R M ) = ^^-e-^. (6) 

The early-time consideration allows us only to determine the region of that has the ridge 
without specifying the qj> distribution of the partons at that time. The partons in the bulk 
and the ridge can take their time to thermalize without affecting in results in Eqs. (jlj) and 
([6]). No rapid thermalization is required. 
From Eqs. (j3J) and ([6]) we obtain 

R{p T A) = R(ptM4>), 

where 

R(Pt) = ^~ VTlT ' (l - e~ PT/T ") , (7) 

0(0) = 0($-|0|) + e($-|7r-0|), (8) 
111 AT 

— = -- — = — AT = T' - T (9) 

rp„ rp rp, rprp,1 L ^ ± WJ 

Experimental study of the ridges in triggered events in central Au+Au collisions shows that 
AT is in the range of 40-50 MeV 0. We take AT = 45 MeV in the following. 

We note that the values of C in Eqs. ([3]) and (jSJ) are the same so that R{pr) ^ as 
Pt — > 0. Since partons at qx = do not flow, only the q T > partons can lead to (f> 
anisotropy, so B + R cannot be distinguished from B at pt = 0. There being no data on the 
ridge at px ~ 0, setting C to be the same in Eqs. ([3]) and (J3J) is a reasonable assumption. 



3 Elliptic Flow 



The second harmonic in the d> distribution is [10 



J 27r dcj) cos 2(f) dN/p T dpTC 
Jq W d(p dN/pTdprdcj) 



Mpt) = (cos20) = 2n - ;Ar; ; . (10) 



When we substitute for dN / pTdpxdcf) the general expression in terms of B(pr) arid R(pr, 
we obtain 

fdcp cos 2(j)R(p T ,(j)) 

V2{PT,0) ^ 



Jd(f)[B(p T ) + R(p T , 
sin 2$ (b) 



7rB(p T )/R(p T ) + 2<$>(b) 

When Eq. (J7|) is rewritten as 



R(Pt) = %-^~ VTlT (f TlT " - l) , (12) 
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we have 



B(pt)/R(pt) = (e^ T " - (13) 

which is large at small p?, and dominates in the denominator of Eq. (1TT1) . Thus we obtain 
for p T < 0.5 GeV/c 

4(Pt,6)^ ^ sin 2S(6). (14) 

This is a simple formula that describes both the p<r and centrality dependences at small pt- 
Note that R{pr) i n Eq. (I12p is proportional to B{px) in Eq. (j3J), so C 2 as well as the 
exponential factor are canceled out in Eq. ffl3|) . Such factors can depend on centrality, but 
are irrelevant to v^ipr)- T" may have some centrality dependence, but it would be mild 
compared to that of the factor sin 2$ (b) that depends on b strongly and explicitly. The 
peak of sin2$(6) at $ = 7r/4 corresponds to b = \^2Ra, which in turn corresponds to ~50% 
centrality for any colliding nuclei. The data indeed shows a maximum of v 2 at that centrality 
for small px [H] (see Fig. 1). From the 7r + distribution at 40-50% [12], a fit of the exponential 
behavior in the range 0.5 < pr < 2.7 GeV/c gives T = 0.287 GeV. Using that value of T 
and 0.045 GeV for AT, we obtain from Eq. (Q 

T" = 2.12GeV. (15) 

It then follows from Eq. (|I4"1) that at small px the maximum value of (at 50% centrality), 
divided by pr, is 

V ^ h) = -L = 0.15 (GeV/c)-\ (16) 

p T txT" y 1 ' ' v ; 

which agrees with the data that shows v-i = 0.075 at pt = 0.5 GeV/c. This is the first 
quantitative demonstration that v-i and ridge phenomenologies are related. For higher pt 
the approximation given in Eq. (fl4"l) is not good enough, so the full expression in Eq. (TTBl 
should be used in fTTTj) . The result is shown by the thick line in Fig. 1 (left) for 50% centrality, 
and agrees well with the data [IT] indicated by the blue stars for 40-50% centrality. 

It should be remarked that data on the ridges are for p? > 2.5 GeV/c [2]; in that pt 
range the distribution can be described by either Eq. (jSJ) or ([7j), giving essentially the same 
T' with errors large enough to render AT =40-50 MeV. Our use of AT =45 MeV for smaller 
Pt in the absence of relevant data is an optimistic extrapolation that allows us to make an 
immediate connection between the ridge data and v 2 . 

For proton production the same formalism of recombination applies, as for pion, and the 
ridge consideration is the same. To take the mass effect into account we do what is usually 
done for hadronization, i.e., replacing all exponential forms, such as e~ PT ^ T for the bulk, by 
exp[— (mj — m)/T], where tut = [p\ + m 2 ) 1 ^ 2 and m the proton mass. It then follows that 
at small Pt the v 2 formula for proton is 

^(Rr,6) = 2^siii2$(6). (17) 
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Figure 1: Calculated v^ipr) and v^ipr) for 50% centrality are indicated by the thick blue 
lines. The data are from [TT] for Au-Au collisions at 200 GeV, together with the light lines 
from that reference. 



With the same T" as in Eq. ( 1T4"1) we have universal behavior of the slope of v^(pt, b) versus 
the transverse kinetic energy, Et, which is pt for pion and p T /2m for proton 



dv% (pr, b) 



3En 



p T <0.5 



7rT 



-sin 2 $(6) 



that is independent of the hadron type h. This universality is a consequence of T" being 
a property of the partons in the ridges before hadronization. Thus it trivially satisfies the 
quark number scaling at low pt, for which both v\ and Et are scaled by the quark number 

/>„ mm. 

For < pt < 1.2 GeV/c we show by the thick blue line in Fig. 1 (right) the calculated 
v 2{Pt) for $ = vr/4 corresponding to 50% centrality; it agrees well with the data for 40-50% 
[TT] . For more peripheral collisions Eq. (|17p is inadequate to describe proton production be- 
cause less thermal partons render the joint wu<i-distribution not factorizable and the effective 
value of T in the bulk is reduced. That leads to a decrease of B{j>t) / R(pt) and the second 
term in the denominator of Eq. (iTTj) cannot be ignored. The consequence is that v^ipr) 
continues to increase for larger b, a feature that will be investigated in detail in the future. 
For the main range < b < V2R A our simple treatment described above is adequate. 

If we define the scaled impact parameter to be (3 = b/2R^ the explicit factor that 
describes the centrality dependence is 



f(P) = sin2(cos- i /5) 



(19) 



which is plotted in Fig. 2. It agrees with the Au-Au data on v% [13], normalized to unity at 
maximum. Since /(/3) does not depend on Ra explicitly, one expects it to agree with the 
Cu-Cu data also. Evidence is given in Ref. 
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The next question concerns the dependence of T" on energy and density. Properties of 
the ridge have so far been studied with the help of triggers, and our knowledge about ridges 
with unspecified semi-hard parton momenta is rudimentary. Our assumption that T" is a 
constant is a reasonable first approximation; however, we expect it to be improved upon 
further investigation. 

The above consideration applies only at midrapidity. In forward region the semi-hard 
scattering involves partons at larger x, but not too large because of recombination |15j . 
Lower parton density means that there are less phantom jets, so the ridge effect is reduced 
with increasing 7], resulting in diminishing t>2. 



4 Summary 

We have shown that the observed azimuthal anisotropy at low can be reproduced by a 
consideration of the ridges whose formation is due to semi-hard scattering of partons. The 
hadrons in the ridges are produced by recombination of enhanced thermal partons near the 
surface of the medium. The centrality dependence is described mainly by a geometrical 
factor, sin 2$ (b). The initial slope in the dependence of t>2 is essentially 1/ttT" for both 
pion and proton, where T" is related to parameters in ridge phenomenology. This is the 
first time that a connection is made between properties of elliptic flow and of ridges due to 
semi-hard scattering. No fast thermalization is required, and no hydrodynamical flow from 
early times is assumed. We do not exclude the relevance of hydrodynamics at some later 
stage of the expansion, but it is not needed to derive the main features of i> 2 from ridge data. 

This work was supported, in part, by the U. S. Department of Energy under Grant No. 
DE-FG02-96ER40972. 
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